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<N ■ ABSTRACT 

<^ . We present a simple numerical model of the plasma flow within the open field 

', line tube in the pulsar magnetosphere. We study how the plasma screens the 

rotationally induced electric field and maintains the electric current demanded 
by the global structure of the magnetosphere. We show that even though bulk of 
the plasma moves outwards with relativistic velocities, a small fraction of particles 
is continuously redirected back forming reverse plasma flows. The density and 
composition (positrons or electrons, or both) of these reverse flows are determined 
q by the distribution of the Goldreich- Julian charge density along the tube and by 

£3 I the global magnetospheric current. These reverse flows could significantly affect 

the process of the pair plasma production in the polar cap accelerator. Our 
simulations also show that formation of the reverse flows is accompanied by the 
generation of long wavelength plasma oscillations, which could be converted, via 
\0 [ the induced scattering on the bulk plasma flow, into the observed radio emission. 

^ . Subject headings: plasmas-(stars:) pulsars: general 

o : 

O ■ 1. Introduction 

> : 

The pulsar activity is believed to be associated with the generation of relativistic 
electron-positron plasma near the magnetic polar caps. This plasma flows along the open 
field line tube and eventually escapes from the magnetosphere forming a relativistic pulsar 
wind. Well within the light cylinder, the plasma currents do not distort significantly star's 
magnetic field therefore in the frame corotating with the star, the plasma just moves along 
the axis of the rotating static dipole. An important point is that this motion could by no 
means be considered as free streaming because the basic electrodynamics dictates the charge 
and current densities at each point of the flow. 

First of all the charge density in the plasma should be equal to the local Goldreich- Julian 
charge density 
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where B is the local magnetic field, VL the angular velocity of the neutron star. This condition 
is a generalization of the standard condition of the quasi-neutrality: deviation of the local 
charge density from pgj results in a longitudinal electric field, which redistributes the charges 
to establish p = p G j. The necessary field is weak in the sense that the corresponding potential 
is small as compared with the total rotationally induced potential; this is because the energy 
of the secondary plasma particles is small as compared with the total potential. 



Scharlemannl (119741 ) and I Cheng fc Rudermanl (119771 ) noted that in the plasma flow along 



the curved magnetic field lines, some difference should be maintained between the velocities 
of the electrons and positrons. The condition p = pgj implies n + — n~ oc B cos 9, where n^ 1 
are the number density of positrons and electrons, correspondingly, 9 the angle between the 
pulsar rotation axis and the magnetic field. Continuity of the particle flow implies n^V^- oc 
B, where V ± are the average velocities of the positrons and electrons, correspondingly. One 
sees that the average particle velocities should vary along a curved magnetic field line in 
accordance with the variation of 9. In a highly relativistic flow, are close to c therefore 
even a small variation of 9 implies a significant variation of the Lorentz factor of the positrons 
and electrons. Since the energy spread of the secondary particles is large, the electric field, 
which ensures the adjustment of the charge density, easily shifts the low-energy tail of the 
distr ibution function to the negative-m omentum domain thereby forming a return particle 
flux (jLyubarskijlll992l ; lLyubarskiilll993al ) . Therefore one can expect counterstreaming plasma 
flows in the open field line tube. 

If the charge density in the flow is adjusted via sending "extra" charges backwards, an 
electric current appears in the flow. However, the electric current is not a free parameter, 
which could be adjusted to maintain the necessary distribution of local charges. The current 
in the open field line tube is dictated by the global structure of the magnetosphere. The 
magnetic field lines are bent backwards with respect to the rotation direction in order to allow 
the plasma to stream along them at a velocity smaller than the speed of light even though 
the rotation velocity becomes superluminal beyond the light cylinder. The current should 
be distributed such that the necessary magnetic field is maintained (formally, the current 
along any magnetic field line is deter mined from the condition of the smooth transition 
of th e flow through the light cylinder (jlngrahaml 1 19731 ; IContopoulos et al.l Il999l ; iTimokhin 
20061 )). Generally the current required by the global structure of the magnetosphere is not 
matched with the current established in the course of the charge density adjustment. This 
would result in some induction electric field, which ensues additional redistribution of the 
charge flows. These considerations show that the self-consistent distribution of charges and 
currents in the open field line tube could be maintained only in a non-trivial, generally 
unsteady, multi-stream state. 
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As the first step to understanding of the electric charge and current adjustment in 
the pulsar magnetospheres, we study here a simple model of the plasma motion within a 
narrow tube with a distributed background charge imitating the Goldreich- Julian charge. 
The plasma production is not addressed here; we just assume that the relativistic electron- 
positron plasma with large enough density (significantly larger than the Goldreich- Julian 
density) is injected continuously into the tube. The aim of this study is to keep track of how 
the plasma screens out the background charge and at the same time maintains the current 
imposed at the outer edge of the tube. The system evolution is studied numerically by the 
particle-in-the-cell (PIC) method. 

It will be shown that the electric charge and the current densities are adjusted in the 
plasma flow such that the condition p = pqj is fulfilled at any point and at the same time 
the current imposed by the boundary conditions is maintained. This justifies the MHD 
approach in studying the global structure of the pulsar magnet osphere. The charge and 
current adjustment is achieved via formation of relatively weak backward flows. This could 
have important implications for the models of the plasma production in pulsars because the 
plasma flow from the magnetosphere onto the polar cap could alter significantly the electric 
structure of the polar gap accelerator. Moreover, long wave plasma oscillations are excited 
in the flow in the process of continuous redistribution of charges. The induced scattering 
of these oscillations by the bulk plasma flow could give rise to the observed pulsar radio 
emission. 

The paper is organized as follows. In the next section, the numerical method is outlined. 
In section 3, we present the results of PIC simulations of the relativistic plasma flow within 
a narrow tube with a distributed background charge. In section 4, we discuss the plasma 
oscillations observed in the simulations. In sect. 5, we generalize the model adding an 
arbitrary current at the outer edge of the tube. The conclusions are given in section 6. 



2. Quasi- ID electrodynamics 



In this section, a simple generalization of the one-dimensional PIC code is described, 
which pinpoints the basic qualitative properties of the plasma motion within the open field 
line tube. Since the particles move along the magnetic field lines like beads on a wire, their 
motion is governed only by the longitudinal electric field, 



— = ±eE z 
dt 



(2) 



Ther efore one can use the parti cle mover from the electrostatic one-dimensional PIC code 



[e.g., 



Birdsall fc Langdonl (Il99ll )). The problem is that in our case, the electric field is not 
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one-dimensional because the plasma moves within a narrow tube; therefore one has to solve 
Maxwell's equations within the tube. We take into account the effects of finite transverse 
dimensions by substituting 1/R, where R is the tube radius, for the transverse derivatives 
in Maxwell's equations. The same substitution is sometimes made in the Poiss on equation 
when studying steady state flows in the open field line tube (IFawley et al.lll977l ). 



Maxwell's equations in the frame rotating with the star could be wr itten as customary 
Maxwell's equations with specific source terms (e.g.. IFawley et al.l (119771 )) 
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Well inside the light cylinder, fir < 1, one can take p = pcj] j = 0. 

Evolution of the fields is governed by Eqs.flH]) and dHj); Eqs. and (jlj) are satisfied 
identically provided the initial conditions satisfy these equations and the charge is conserved. 
Within the cylindrical tube, the evolution equations are written as 



BE* 

dt 
dE r 

~dt 

dB„ 



1 drB v 4tt 

r or c 
dB, 



<p . 



dz 
dE r 



(7) 



+ 



dE x 



dt dz dr 

Here B v is the azimuthal magnetic field created by the currents in the plasma. Well within 
the light cylinder, this field is very weak as compared with the pulsar magnetic field therefore 
it does not affect the particle motion; the particles move in the z direction along the pulsar 
magnetic field lines. We assume that the walls of the tube are highly conductive so that 
Eg = at the boundary of the tube, r = R. Substituting the transverse derivatives by 1/R, 
we get a one- dimensional set of equations 



dE K B„ 



An 



dt 
dE r 

~dt 

dB ir 



R c 



dB v 

dz 
dE r 



E-, 



dt dz R ' 

which is qualitatively equivalent to the original Eqs. (j7]). Here we took into account that 
E g goes to zero at the wall so that dE z /dr =>■ —E z /R. Eqs. (JHJ) could be reduced to the 
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one-dimensional, inhomogeneous Klein-Gordon equation. At the scale much less than the 
radius of the tube, these equations are reduced to the one-dimensional Maxwell equations. 
At larger scales, they correctly describe the exponential decreasing of the electric field from 
a localized charge due to the presence of conducting walls. For boundary conditions, one 
sets E r at the bottom edge of the tube and B v at the outer edge. 

We start from the empty tube therefore B^it = 0) = whereas the initial electric 
field is found from the solution to the Gauss equation with p = 0. With the same 
prescription for the transverse derivatives as in equations (|S]), one can write E r (t = 0) = <& / R; 
E z (t = 0) = —dQ/dz; where the potential $ is found from the modified Poisson equation 

With the initial conditions thus found, the evolution equations (jHJ) are solved by the 
leap-frog method: 

At 47T 

E z (z, t + At)- E z (z, t) = -ft B <pi z > t + 0.5At) - —Atj(z, t + 0.5 At) 
E r (z + 0.5Az, t + At) - E r (z + 0.5Az, t) 

= ~ a7 ^ z + Az > t + °- 5At ) ~ B ^ z > t + °- 5At )] ; ( 10 ) 

B v (z, t + 0.5At) - B v (z, t - 0.5At) 

= --r- - [EJz + 0.5A^, t) -EJz- 0.5Az, t)] - —EJz, t). 

Az it 

The current density, j , is found in each step by calculating the charge transferred by particles 
through each grid face. We use the first order particle weighting, i.e. the particle is assumed 
to be a cloud of uniform density, the cloud width being equal to the grid size. 

Stability of the scheme is analyzed as usual by writing down the dispersion relation for 
the homogeneous system: 

2R\ 2 . 2 fuAt\ 1 (2R\ 2 . 2 fkAz\ 

sm -S- = 1 + sm -S- • ( n ) 



At y v 2 / V A <v V 2 

One sees that under the condition 

< At > 2 * (12) 

the frequency u; is real, which implies stability. In the simulations, we used At = 0.5Az. 



- 6- 



3. Screening of the rotationally induced electric field 

The rotationally induced electric field implies specific source terms in Maxwell's equa- 
tions written in the rotating frame of reference. Deep inside the light cylinder these terms 
are reduced to the Goldreich- Julian charge density of equation (JTJ), which could be consid- 
ered as a background charge distributed within the tube. Our main interest here is with the 
compensation of the background charge in the relativistic plasma flow. In the static system, 
the particles of both charge species would easily be redistributed to provide total neutrality. 
The situation is less trivial when the plasma flows through the tube. We investigate the 
simplest case when the plasma density in the flow, n, is large enough to compensate the 
background charge, en ^> p GJ . 

For the sake of simplicity, we assume that the radius of the tube, R, is constant. In the 
pulsar magnetosphere, the open field line tube expands with the altitude, which results in 
decreasing of the particle density as well as of the charge and current densities. However, the 
Goldreich- Julian charge density ((Tj) decreases by the same factor because the pulsar magnetic 
field decreases in the expanding tube. Therefore only variation of the angle between the 
magnetic field and the rotation axis, but not the tube expansion, results in the discrepancy 
between p and pgj- Thus the plasma flow in a tube of a constant radius with a varying 
background charge seems to be a reasonable model for studying screening of the Goldreich- 
Julian charge in the open field line tube. Note also that only variation of the background 
charge along the tube results in the nontrivial behavior of the plasma. If pgj = const, the 
charge density would be adjusted already at the bottom of the tube and then the acquired 
charge would be just transferred by the flow further into the tube. As we are interested in a 
continuous charge redistribution in the plasma flow along the tube, we take the background 
charge to be varying along the tube. The simplest choice is a linearly growing background 
charge density 

Pgj = poz/k; (13) 

where Iq is the length of the tube. Then the system is neutral at the bottom of the tube 
whereas at higher altitudes, charge neutrality may be maintained only via redistribution of 
charges in order to compensate the growing background charge. 

We took the radius of the tube much smaller than the tube length, R = /o/30, which 
means that different parts of the tube are not connected electrically. The initial electric field 
was found from the solution to the Poisson equation with E z = at both edges of the 
tube. One can see that far from the edges of the tube, z ^> R and Iq — z 3> R, the solution to 
equation (jUJ) with the background charge density ffTBT) is $ = —AnpoR 2 z/lo. Then the initial 
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longitudinal electric field is constant far from the edges of the tube, 

E z (t = 0) = A7ip R 2 /l ; (14) 

whereas the transverse electric field increases linearly from the bottom to the upper edge of 
the tube, 

E r (t = 0) = -A7rp Rz/lo. (15) 

We start injecting plasma from the bottom into the empty tube. At the bottom bound- 
ary, we use the condition E r (z = 0) = 0, which follows from the continuity of the transverse 
component of the electric field. At the upper boundary, we take B^z = l ) = 0. Then 
the plasma moves as within an infinitely long tube until the flow approaches the outer edge 
of the tube by the distance R, i.e. while the plasma is disconnected electrically from the 
outer boundary. Electrons and positrons are injected continuously with the same density 
and distribution function; the last is shown in Fig. 1. The pair density in the flow is n = 10 
pairs per cell. We choose po — O.len so that the screening of the background charge would 
require redistribution of a small fraction of the particles. The electron charge, e, was chosen 
such that uj p Az/c = 0.3 whereas the radius of the tube was R = lOOOAz. Introducing is 
the characteristic plasma frequency, u p = y^87ie 2 n/m (recall that n is the number density 
of pairs therefore the total particle number density is 2n), one sees that u p R/c = 300. With 
the chosen parameters, the total potential drop in the empty tube is 

e$ = AnepoR 2 = —muj 2 p R 2 = 4500mc 2 ; (16) 

which is much larger than the energy of the injected particles. Therefore the first injected 
particles are accelerated to high energies. However within the plasma, the electric field is 
nearly shorted out. Compensation of the background charge is provided by decelerating and 
sending back "extra" charges (electrons in our case) therefore within the plasma flow, a small 
electric field is maintained sufficient to redirect the slowest particles. 

Evolution of the fields and the particle phase space distribution is shown in Figs. 2 and 
3. One sees that an electron backflow is indeed formed and the plasma efficiently screens the 
longitudinal electric field. At the bottom of the tube, a relatively large electric field arises 
with the total potential drop roughly 10mc 2 /e. This is because in the injected plasma, there 
are no particles with the energy less than lOmc 2 so that a "double layer" is formed in order to 
shift the low energy tail of the electron distribution to the negative momentum region. In the 
remainder part of the tube, the field is very low so that the total potential drop is comparable 
with that near the bottom edge. The distribution of the electric potential along the tube 
is determined by the distribution of the background charge and by the particle distribution 
function in the low energy part. The background charge dictates how many electrons should 



- 8- 




Fig. 1. — The distribution function of the injected particles (thick line), as well as the 
distribution functions of positrons (thin line) and electrons (dots) found in the region 0.25/o < 
z < 0.75/q at the end of simulations, at t — 3.25/ /c 
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be removed from the flow and the potential is adjusted such that the necessary amount of 
electrons is redirected back. 

As the plasma moves farther into the tube, the increasing number of electrons is sent 
back in order to compensate the increasing background charge. The increasing reverse flux 
of electrons implies the increasing electric current flowing in the plasma and therefore the 
increasing azimuthal magnetic field. Note that the plasma screens immediately only the lon- 
gitudinal electric field whereas the transverse electric field remains nonzero; it only decreases 
as compared with the initial transverse field. As the transverse field varies along the tube, 
the curl of the electric field arises, which is balanced, according to the induction law, by the 
growing magnetic field. 

In order to check whether a steady state solution is possible, at least in the average, we 
extend the tube to z > l Q assuming that the background charge density remains constant at 
z > l : 



One can anticipate that in the region with p G j — const, the reverse plasma flow is not formed 
therefore this region would not affect the region z < Iq so that a steady state solution would 
be eventually formed there. The simulations show that this is indeed the case (see Figs. 4 
and 5). The transverse field gradually decreases to zero and a steady state is formed with 
E = and p = pgj- m the course of time, the steady zone extends further into the tube. 

One can easily estimate the current, and therefore the azimuthal magnetic field, estab- 
lished in the steady state. In the region z > Iq, the background charge density is constant, 
Pgj — Po] the particle redistribution does not occur there and the reverse particle flux is not 
formed. Then the current may be presented as j = ec(n + — where n + and n_ are the 
number densities of positrons and electrons, correspondingly. Here we assume that all par- 
ticles move with the speed of light. In the pulsar magnetospheres, the particle distribution 
is wide, A7 ~ 7 ~ 100, and in the presence of the longitudinal electric field, the distribution 
function is just shifted in the momentum space as a whole so that the fraction of the particles 
with the modulo velocities significantly less than the speed of light remains negligibly small 
in any case. Taking into account that charge neutrality implies e{n + — nJ) = po, one finds 
finally j = cp^. When the steady state is achieved, the magnetic field approaches 




z < l ; 

Z > l . 



(17) 



B vfi = AnRp . 



(18) 



This estimate assumes that the electric field in the tube is screened completely. However, 
some small field remains in order to decelerate the necessary amount of electrons and send 
them backwards. Therefore the charge density in the plasma is a bit less than the background 
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Fig. 2. — The phase space of electrons (a), positrons (b) and the distribution of the nor- 
malized fields, e z = E z l$/ (Anp R 2 ) (solid), e r = E r /(4irp R) (dotted), b = B^/ (A-Kp R) 
(dashed), at t — 0.25/ /c. The fields are smoothed at a scale 2nc/u p . 




Fig. 3. — The same as in Fig. 2 at t = 0.67/q/c. 




Fig. 4. — The same as in Fig. 2 at t — 2.4/ /c. 
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Fig. 5. — The same as in Fig. 2 at t = 3.6/q/c. The steady state is formed at z < l 
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charge density, which implies that the magnetic field remains a bit less than that of equation 
flUD; see Fig. 5. 



4. Plasma oscillations excited in the course of the charge adjustment 

Inspection of Figs. 4 and 5 shows that strong fluctuations of the longitudinal electric 
field are present in the region < z < l where the background charge density varies. In Fig. 
6, we show the power density spectrum of the fluctuations, P = (87i)~ 1 \E Zy k\ 2 , in this region 
as well as in the region z > Iq. One sees that in the last region, where the background charge 
density is constant so that additional redistribution of charges does not occur, only short- 
wavelength fluctuations are present, which could be attributed to the numerical noise. In 
the region < z < lo, the fluctuation spectrum is extended to the long- wavelength domain, 
which means that the charge redistribution is accompanied by generation of long-wavelength 
plasma oscillations. 

The longitudinal fluctuati ons in the relativi stic,one-dimensional plasma are described 



by the dispersion relation (e.g. iLyubarskiil (119961 )) 



2 f f+ip) + f-(p) i 1 / 1f rt 

U >J 7 3(^-H 2 P= ; ( 9) 

where f± are the distribution functions of electrons and positrons, correspondingly. The 
characteristic frequency of the long wavelength oscillations, 

( f (U + f-)d P V /2 

p u — ~ 3 — ) ' ( ' 

is determined by the low energy particles. In Fig. 7, we show the numerical solution to the 
dispersion equation (fT9"|) for the distribution functions found in our simulations (and shown 
in Fig. 1); one sees that u pa 0.06u;p in this case. According to Fig. 6a, the maximum 
of the fluctuation energy, kP, occurs at k ~ 0.05u p /c, which is close to ujq/c. In order to 
check this scaling, we performed the same simulations with uj p three times larger than in the 
simulations shown in Figs. (2-5). In this case, the maximum of the power spectrum shifted 
three times towards a lower k, as one should expect. 

Origin of these fluctuations is not very clear. For steady state boundary conditions, 
one can easily construct a completely time-independent state. The amount of particles redi- 
rected at each point is determined by the spatial variation of the Goldreich- Julian density. If 
the injection is time-independent, the potential necessary to redirect the required amount of 
particles is also time-independent. In principle, the obtained state is unstable with respect 
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to two-stream instability because the momentum distribution of the injected particles grows 
with p at small p > whereas the distribution of the redirected (and therefore having a neg- 
ative p) particles decreases with p (jLyubarskii 1993a ). However, hardly ever this instability 



operates in the presented simulations because numerical fluctuations of the electric field are 
relatively strong so that the particles are efficiently scattering smearing out any peculiari- 
ties on the distribution function. By this reason, no signs of the two-humped distribution 
function was seen in the simulations. 

Fluctuations presumably arise because we inject equal amount of electrons and positrons 
whereas in the true steady state, the injected flux of electrons should be less by (l/2)c/?o/e 
than the flux of positrons. This is because an electron redirected at a point z = Z\ contributes 
twice into the charge density at z < Z\ so that the total amount of injected electrons should 
be less than that of positrons by amount of the electrons redirected along the whole tube. If 
this condition is not satisfied, extra charges (electrons in our case) are expelled from the flow 
at the scale of the Debye length, which is determined by the frequency ujq for the injected 
plasma. Therefore a "double layer" arises at the bottom of the tube; it is clearly seen in 
Fig. 3 and less clearly (because of small scale) in Figs. 4 and 5. Such a double layer is 
unsteady; it oscillates at the time scale of 1/uJq thus slightly modulating the injected plasma 
flow. This could give rise to the electric field fluctuations in the tube because if the plasma 
flow is somehow modulated, the decelerating electric field should fluctuate in order to keep 
the number of redirected particles fixed (the last is determined by the distribution of the 
Goldreich- Julian charge density, which does not vary). On the other hand, one has to stress 
that the modulation occurs at the frequency ujq for the injected plasma, which is less than the 
frequency ujq for the plasma in the tube because the lowest energy particles in the injected 
plasma have energy 7 ~ 10 whereas in the tube, the particle distribution is shifted to 7 ~ 1. 
So the presented interpretation is not straightforward. 

In real pulsars, the characteristics of the injected plasma depend on the conditions in the 
polar cap region therefore there is no reason to expect that the polar cap accelerator injects 
the particles in such a way that the difference in the positron and electron fluxes is adjusted 
to the structure of the magnetosphere in the light cylinder region (recall that po is referred to 
the Goldreich- Julian density far from the polar cap, somewhere at fir ~ 1). Moreover, the 
plasma production process is by no means steady therefore one can expect that the plasma 
flow in the open field line tube is strongly modulated. Any modulation of the plasma flow 
would result in fluctuations of the decelerating field along the tube thus giving rise to long 
weavelength plasma oscillations in the flow. Therefore we believe that such fluctuations are 
generic in pulsar magnetospheres. Induced scattering of these oscillations o n the bulk plasma 



flow could g ive rise to high brightness temperature pulsar radio emission (ILyubarskiil 11992 



1993U 11990 ) 
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5. Adjustment of the electric current in the plasma flow 

One sees that when the system relaxes to the steady, on the average, state, the current 
(and therefore the azimuthal magnetic field) is established, which is determined by the 
distribution of the Goldreich- Julian charge density. However, there is no reason to expect 
that this current is consistent with the current necessary to maintain the global structure 
of the magnetosphere. As the condition E • B = is fulfilled within the plasma flow, the 
structure of the pulsar magnetosphere satisfies the ideal MHD equations and as such do 
not admit additional constraints on the current because the current distribution in any 
MHD system is uniquely determined. Specifically in the pulsar magnetospheres, the current 
distribution is determined from the condition of the smooth plasma tra nsition through the 



light cylinder (Ilngrahamlll973l ; IContopoulos et al.lll999l ; lTimokhinll2006l ). Let us check what 



happens if the current "demanded" by the global MHD solution differs from the current 
established in the course of the charge density adjustment, i.e., if the field ffTB"]) does not fit 
the MHD structure of the magnetic field. For this purpose, we simulated plasma flow within 
a tube of the length l with a (generally varying) magnetic field imposed at the upper edge of 
the tube. As an initial condition, we used the distributions of the fields and particles found 
at the end of the previous simulations when the steady state was achieved in the region 
< z < l . We imposed the condition of free particle escape from the upper edge of the tube 
and no particle inflow through this edge. Conditions at the bottom edge remained the same 
as in the previous simulations, i.e. the plasma is continuously injected into the tube and 
E r (z = 0) = 0. When we impose the condition B 9 {z = Iq) = B Vj q at the upper edge of the 
tube, nothing changes and the initial state persists forever. However, when we slowly vary 
B^z = Iq) at the upper boundary from B v to another value, -B^i, a discrepancy between 
the current in the plasma and the imposed B v results in the induction electric field, which 
eventually establishes the necessary current. 

In the simulations, the magnetic field at the upper boundary was changed according to 
the formula 

B v {z = l ) = B vA + (B vfl - B Vtl ) exp (-t/r) ; (21) 

In Figs. 8 and 9, we show how the system is relaxed to the state with the azimuthal magnetic 
field B 9t i = 0.5B Vt o. One sees that the electric field arises near the upper edge of the tube 
so that positrons are decelerated until some of them are redirected back thus compensating 
the initial current. Eventually a steady state is formed with the new magnetic field and with 
both electron and positron backflows. The electrons are redirected along the whole tube 
compensating variation of the background charge density. The positron backflow is formed 
near the upper edge of the tube. 



Relaxation to the state with the current larger than the initial one is achieved via 
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formation of an enhanced electron backflow, see Fig. 10. In Fig. 11, we show the state 
with the current having the sign opposite to that of the Goldreich- Julian charge. Such a 
state is achieved by sending more positrons back. So one concludes that the system could 
easily maintain any current imposed at the outer edge, the adjustment being achieved via 
formation of an appropriate particle backflow. 



6. Conclusions 

We have demonstrated that the relativistic plasma flow in the open field line tube of the 
pulsar easily screens the rotationally induced longitudinal (along the magnetic field) electric 
field and maintains any current demanded by the boundary conditions. This justifies the 
MHD approach to the structure of the pulsar magnetosphere provided the plasma density in 
the flow exceeds the Goldreich- Julian density. An important point is that the adjustment of 
the electric charge and the current is achieved via formation of a particle backflow with the 
density roughly pcj/e- Depending on the magnetospheric current and on the distribution of 
the Goldreich- Julian charge density along the open field line tube, the backflow could consist 
of electrons, positrons or both. 

In these simulations, the pair plasma with the density significantly larger than pGj/e 
was continuously injected from the bottom of the tube so that a steady state was eventually 
achieved. In the pulsar magnetosphere, the plasma is assumed to be produced at the bottom 
of the open field line t ube w here primary particles are accelerated in the potential gap (e.g., 



Hibschman fc Arond (1200 11 )). The present study shows that one should expect a particle 
flux coming from above to the gap. When the plasma fills the gap from above, the electric 
field is shorted out, the particle acceleration is terminated and the pair production ceases. 
Presumably this means that the polar accelerator could not operate in the steady state 
regime. One can expect some sort of spasmodic plasma production such that the required 
current is maintained only on average. Future models of the polar cap accelerator should 
take into account the plasma flow onto the gap from the magnetosphere. 

Another implication of this study is that the charge and current adjustment is ac- 
companied by continuous generation of long-wavelength plasma oscillations. As a result 
of the induced scattering of these oscillations in the bulk plasma flow, these oscillations 
could be converted into the electromagnetic waves with the frequency u ~ w PV /7max, where 
1 00 is the Lorentz factor corresponding to the maximum in the particle energy 



spectrum (I L y ub arskiil 1 1 9 9 61 ) . The energy of these waves exceeds the energy of the initial os- 



cillations ~ w PV /7max/ times. One can sp eculate that this p rocess gives rise to the observed 



pulsar radio emission (for more details see iLyubarskiil (119961 )). 



-18- 



I thank Saveli Rabinovich for help in designing the code. Useful discussions with Jon 
Arons and Anatoly Spitkovsky are gratefully acknowledged. This work was supported by 
the Israeli Science Foundation and by the German-Israeli Foundation for Scientific Research 
and Development. 

REFERENCES 

Birdsall, C. K., & Langdon, A. B. 1991, Plasma physics via computer simulations (IOP 
Publishing Ltd, 1991, 479 p.) 

Cheng, A. F., & Ruderman, M. A. 1977, ApJ, 212, 800 

Contopoulos, I., Kazanas, D., & Fendt, C. 1999, ApJ, 511, 351 

Fawley, W. M., Arons, J., & Scharlemann, E. T. 1977, ApJ, 217, 227 

Hibschman, J. A., & Arons, J. 2001, ApJ, 560, 871 

Ingraham, R. L. 1973, ApJ, 186, 625 

Lyubarskii, Y. E. 1992, A&A, 265, L33 

— . 1993a, Astronomy Letters, 19, 14 

Lyubarskii, Y. E. 1993b, Astronomy Letters, 19, 208 

Lyubarskii, Y. E. 1996, A&A, 308, 809 

Lyubarskij, Y. E. 1992, A&A, 261, 544 

Scharlemann, E. T. 1974, ApJ, 193, 217 

Timokhin, A. N. 2006, MNRAS, 368, 1055 



This preprint was prepared with the AAS IATgX macros v5.2. 



-19- 



4 

4, . 1 1 1 




lg(kc/co p ) lg(*c/w p ) 

Fig. 6. — The power spectrum of the fluctuations of the longitudinal electric field in the 
region 0.25/ < z < 0.75/ (a) and 1.5/ < z < 2l (b). 
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Fig. 7. — The spectrum of the longitudinal oscillations in the plasma with the distribution 
functions of electrons and positrons shown in Fig. 1 by the thin solid and dotted lines, 
correspondingly. 
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Fig. 8. — Beginning of the current adjustment. Shown are the phase spaces of the electrons 
(a) and positrons (b) and the distribution of the normalized fields at t = lo/c after the 
magnetic field at the upper (right) boundary of the tube begun to vary according to Eq. 
(J2TJ) with By i = 0.5-B^ o and r = 5l /c. One sees that the decelerating electric field is 
developed near the upper edge and a weak positron backflow is formed; the magnetic field 
is decreased a bit at this stage. 




Fig. 9. — The same as in Fig. 8 at t — 8/o/c. The necessary positron backflow is already 
formed and the magnetic field is nearly two times less than the initial one. 
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